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Are z^5 QSOs found in the most massive high redshift 
over-densities? 
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ABSTRACT 

Luminous high-redshift QSOs are thought to exist within the most m.assive dark mat- 
ter haloes in the young Universe. As a consequence they are hkely to be markers 
for biased, over-dense regions where early galaxies cluster, regions that eventually 
grow into the groups and clusters seen in the lower-redshift universe. In this paper 
we explore the clustering of galaxies around z '^ 5 QSOs as traced by Lyman break 
Galaxies (LBGs). We target the fields of three QSOs using the same optical imag- 
ing and spectroscopy tech niques used in the ESO Remote Galaxy Survey, (ERGS, 
iDouglas et al.l l2009l |2010D , which was successful in identifying individual clustered 
structures of LBGs. We use the statistics of the redshift clustering in ERGS to show 
that two of the three fields show significant clustering of LBGs at the QSO redshifts. 
Neither of these fields is obviously over-dense in LBGs from the imaging alone; a possi- 
ble reason why previous imaging-only studies of high redshift QSO environments have 
given ambiguous results. This result shows that luminous QSOs at z ~ 5 are typically 
found in over-dense regions. The richest QSO field contains at least nine spectroscopi- 
cally confirmed objects at the same redshift including the QSO itself, seven LBGs and 
a second fainter QSO. While this is a very strong observational signal of clustering at 
z ~ 5, it is of similar strength to that seen in two structures identified in the 'blank 
sky' ERGS fields. This indicates that, while over-dense, the QSO environments are 
not more extreme than other structures that can be identified at these redshifts. The 
three richest structures discovered in this work and in ERGS have properties consis- 
tent with that expected for proto-clusters and likely represent the early stages in the 
build-up of massive current-day groups and clusters. 



Key words: galaxies: high redshift 
scale structure of Universe 



galaxies: clusters: general - cosmology: large- 



1 INTRODUCTION 

n 

Our theoretical understanding of the evolution of struc- 
ture appears to be well advanced and has been thor- 
oughly explored thr ough large numerical simulations (e.g. 
ISpringel et al.l 120051 ). Such simulations, using the current 
cosmological parameters determined from observations of 
distant supernovae, the microwave background and other 



Based on observations made at the European Southern Obser- 
vatory Very Large Telescope, Paranal, Chile (ESO programme 
numbers 082.A-0354{A), 083.A-0618(A), 084.A-0251(B), 085.A- 
0444(B) & 087.A-0589(A)(B)). 



tracers, can largely predict the structures that we see in ob- 
servational studies of the relatively nearby universe. 

By contrast, there are currently few observations that 
directly detail the early growth of structures such as groups 
and clusters in the first few Gyrs. Massive clusters can be di- 
rectly observed out t o z '--^ 2 from wide-a r ea optical/IR and 
X-ray surveys (e.g. Bremer et al. 20061: iFassbender et al 



2011 



20121 : 



ay surveys (e.g. itaremer et ai.i izuua: ir'asspenaer et at 
it ISantos et al.ll201ll: ISpitler et al.ll2012l: [Stanford et al 



Zeimann et al.ll2012l ) and proto-clusters to somewhat 



high er redshifts by targeting either powerful r adio sources 
(e.g. iMilev et alll20C)4: Venemans et alll2007l : iHatch et all 
201 11: iKuiper et al.l 1201 li . |2012| ) or submm sources (e.g. 
iCapak et al.ll201ll ). At the highest redshifts we have detected 
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very few systems that we can be confident are the progeni- 
tors of the m assive structures that exist in the curr ent-day 
Universe (e.g. lOta et al.ll2008l : Flb"shikawa et all2012h . There 
is a need to identify and study these individual systems back 
to the earhest stages in their evolution, in order to under- 
stand the interplay between the formation and evolution of 
these structures and the galaxies that they contain. 

While it is possible to identify clustered regions at the 
highest redshifts through spectroscopy of randomly-selected 
deep fields, few have been identified in this manner at 
z > A [the only kn own z > 4 structures are the two in 
iDouglas etla] (120101. hereafter DIP) at z ~ 5, the z ~ 5.7 
structures in Ouchi et ajj (|2005l ) and the z ~ 6 structure in 
[Toshikawa ct al. (2012)]. Carrying out these 'blind' studies 
is potentially observationally expensive, given the volumes 
and therefore number of objects which need to be targeted. 
Clearly, it is beneficial to have more effective methods for 
identifying potential high redshift clustering prior to large- 
scale observational programs. 

It has long been thought that luminous high-redshift 
QSOs should be found in the most massive hosts and haloes 
(|TurneiJ Il99ll ) and that these systems should be associ- 
ated with ove r-densities in the early matter distribution 
iJEfstathiou fc Rees igSS ). Subsequent analytical and com- 
putational modelling indicates powerful high-redshift QSOs, 
hosted in massive dark matter haloes, may trace some of 
the most over-dense regions in the early Universe, and 
that these regions are li kely to be the progenitors of mas- 
sive low-redshift clusters JSpringel et al.ll2005l:lLi et al.l2007l: 



Trenti fc Stiavellill2007l:lDi Matteo et al.ll2008l:lAngulo et all 



20121 '). The most recent of these works predict that the 
current-day descendants of structures hosting luminous 
high-redshift QSOs have a range of masses from ~ 6 x 
10^'^ Mq to more than 10^^ M© (i.e. the most massive low-z 
clusters), with the typical mass approaching 10^* Mq. 

Prompted in part by these predictions, several ob- 
servational stud i es (Istiavelli et al.||2005 : Zheng et al.l 120061 : 
iKim et al.ll2009l : lUtsumi et al.ll2010l l have probed the envi- 
ronments of z Js 5 QSOs for galaxies at the same redshift, in 
order to identify present day clusters in the earliest stages of 
formation. These studies have resulted in an ambiguous pic- 
ture, with four of the QSO fields studied showing evidence 
for clustering and three showing no evidence for significant 
over-densities. In part this is due to the methods used. Many 
of the studies have targeted z ^ 6 QSOs and then searched 
for over-densities in the number counts of dropout or Ly- 
man break galaxies (LBGs), using photometric data with 
little or no spectroscopy. While the LBG technique can, and 
does, select galaxies out to the highest redshifts, due to the 
increasing luminosity distance, brighter ground-based sky 
against which to selec t and the rapidly-evolving luminosity 
funct ion with redshift (jBouwens et al.ll201ll : lLorenzoni et al.l 
[2013), it is significantly more challenging to work at z ~ 6 
than at z ~ 5 - despite the difference in look-back time be- 
ing just ~ 240 Myr. Civen the variation in number counts 
at faint levels and on arcminute scales, searching for an ex- 
cess of sources in the near-environment of a QSO is likely 
to be a blunt test of clustering of galaxies at the QSO red- 
shift. The broad-band drop-out selection typically identifies 
sources over a range of redshifts ( Az ~ 0.5 — 1) and any clus- 
tering is likely to be on scales of Az '^ 0.1 or less. Thus, an 
over-density would have to be sufficiently strong to be reli- 



ably detected as a perturbation in the galaxy number counts. 
Given that the surface density of z ~ 5 — 6 drop-outs seen 
in d eep ground-based data i s typically 1 per few arcmin'^ 
fe.g. lLehnert fc Bremeiil2003r ). an excess of just a handful of 
objects within Az ~ 0.1 of a QSO's redshift is a significant 
signature of clustering. However, this signal would be lost 
in the background number counts without comprehensive 
spectroscopic follow up. 

In this work we target the fields of three z ~ 5 QSOs, us- 
ing an identical approach to tha t which we used in the ESO 
Remote Galaxy Survey (ERGS: lOouglas et al.|[2009 . 2010 ). 
In ERGS we spectroscopically confirmed 70 z ~ 5 LBGs 
in 10 ^ 45 a^cmin ^ non-QSO fields using VLT/F0RS2 
(JDouglas et al.ll2009l . DIO). While the majority of the ERGS 
fields showed no evidence for redshift clustering (typically 
there were zero or one LBG in each Az = 0.1 bin between 
4.6 < z < 5.6), two fields contained 'spikes' in their red- 
shift distribution indicating significant clustering of galaxies 
into structures. The multi-wavelength properties of those 
galaxies and structures have been probed in several p apers 
(|Stanwav et al.ll2008bl . [201ol : iDavies et al.ll2010l . 12013 ). The 
purpose of this paper is to make a fair and robust compari- 
son between fields which are either centered on a luminous 
QSO or not. By using a similar observing strategy for the 
QSO fields as in ERGS, not only can we be sure that if these 
QSOs are found at the centre of over-densities of LBGs then 
we will identify this large scale structure, but we can also 
compare the strength and character of any clustering to that 
seen in the 'blank' sky fields of ERGS. 

All opti cal magnitudes in th is paper are quoted in the 
AB system (|Oke fc GunnI 119831 ) and a ACDM cosmology 
with Ho=71 kms"^Mpc"\ f2M=0.3 and f2A=0.7 is used 
throughout. 



2 OBJECT SELECTION AND OBSERVATIONS 

The three QSOs targeted in this study were chosen from the 
~ 20 SDSS QSOs with Decl. < 10° and spectroscopically- 
confirmed redshifts of between 4.9 < z < 5.2 prior to 2009. 
The three objects were selected to give an available target 
at any time during the year, i.e. to be ~ 8 hours apart in 
RA, and to have a small range in UV luminosity of no more 
than an order of magnitude. The redshift range was chosen 
to match the most sensitive region of the ERGS selection 
function (see fig. 7 of DIO). All three of the quasars chosen 
are very luminous with absolute magnitudes of between -28.1 
to -29.5. One of these fields (J0338-I-0021) has a second, less- 
luminous, non-SDSS QSO at a sim ilar redshift in the same 
FORS2 field of view (|Diorgovski et al. 2003). However, this 
was not taken into account when selecting this field. 

Three ~ 7' x 7' regions centered on the QSOs were 
imaged in R, I and z using FORS2/VLT and subsequently 
reduced in an analogous manner to the ERGS fields (see 
Table [1} ■ Sextractor (Berlin & Arnouts 199Q) was used in 
dual mode to create 7-band selected catalogues of objects in 
each of the fields, z ~ 5 LBG candidates were then selected 
using the same techniques and colour selection criteria as 
those used in DIO, but with minor alterations intended to 
increase the efficiency of targeting galaxies at the QSO red- 
shift. A primary cut oi R ~ I > 1.3 was applied to select 
LBG candidates as, in the earlier work, this proved to be 
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successful at identifying z ~ 5 galaxies (fig. 7 & 8 of DfO). 
Given the uncertainties in the photometry, some objects at 
4.9 < z < 5.0 may be missed using this cut, especially if 
the object has sig nificant Lyg emission c ontributing to the 
_R— band flux (see IStanwav et al.l l2008bl . for a discussion). 
One of the QSOs, J2130+0026 at z = 4.950, has a colour 
in the FORS2 fihers of _R - J = 0.9 (see Table 1) due to 
Lya+NV emission contributing to its _R-band emission (an 
observed equivalent width of 2020 A). While LBGs do not 
have such large equivalent widths (see fig. 12 in DIG), there 
is some chance that objects in this field, at the redshift of 
the QSO, could be missed by the primary R — I selection. 
Therefore, our colour constraint was relaxed to include bluer 
objects with R — I > 1.0. 

The multi-colour images of each candidate object were 
then examined by eye and, because of our limited number of 
spectroscopic masks, prioritised for potential inclusion fol- 
lowing the prescription outlined in DIO. As in the previ- 
ous work, we avoided using an explicit / — 2 cut. However, 
because we found a trend of increasing I — z colour with 
redshift between 4.6 < z < 5.6 in ERGS (as expected), 
and in this work are specifically searching for galaxies at 
the lower end of this range, we gave higher priority to ob- 
jects with I — z < 0.5. Approximately 20 per cent of the 
spectroscopically-confirmed ERGS galaxies at 4.8 < z < 5.2 
were redder than this. Hence, such systems were not rejected 
outright. The priority assigned to an object was also de- 
pendent whether or not its (ground-based) morphology was 
consistent with those of previously confirmed LBGs (unre- 
solved or barely resolved) and how close its colours were to 
the stellar locus (objects brighter than / ~ 25 and red R~ I 
and I — z colours are almost invariably stars and low red- 
shift galaxies). We also considered whether or not an object's 
photometry was uncompromised by crowding, confusion or 
scattered light from a nearby brighter object. The highest 
priority objects were classified as either priority 1 or 2, the 
first category indicating that the photometry was consistent 
with the source being at z ~ 5 and that there appeared 
to be no issues with the photometric measurements. Those 
objects that had I — z > 0.5 or were brighter than typical 
2: ~ 5 galaxies were assigned as priority 2. 

Two spectroscopic masks were created for each field. As 
with DIO, sources were placed on the masks with the higher 
priority object having preference where dispersed spectra of 
two or more objects contended for the same detector area. 
In general, there appeared to be little obvious 2-dimensional 
spatial clustering of the high priority targets, which would 
compromise the number of these sources placed on the masks 
through slit contention. Each mask was observed for 12800 
seconds split into 20 exposures dithered along the slit, us- 
ing the 3001 grism and the OG590 order sorting filter. Ap- 
proximately sixty objects were targeted in each field - the 
numbers of priority 1 and 2 objects targeted in each case are 
given in Table [21 



3 RESULTS 

Each spectrum was visually inspected by four of the authors 
(KH, MNB, LJMD and ERS) in order to classify the target 
and/or determine its redshift. For those objects identified as 
LBGs, we used the same classification as in DIO, defined as 
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Figure 1. An example LBG spectra with strong Lyman alpha 
emission giving a definite redshift (grade A) . Note that the noisy 
vertical fines in the 2D plot are the residuals from subtracting the 
OH- emission from the sky. The sky fines are plotted underneath 
for comparison. 



follows. Those with a clear continuum break and/or emission 
line indicating a precise redshift were denoted grade A. The 
resolution of each spectrum is ~ 400 kms~^ at z=0, but 
in regions where a continuum break fell in spectral regions 
containing strong sky-line residuals, the break could only be 
located to an accuracy of ~ 1000 kms~^ at z=0. LBGs with 
no emission lines but with breaks falling in such regions were 
denoted grade B. 

As well as those objects spectroscopically confirmed as 
LBGs, the spectroscopy identified a similar range of contam- 
inating objects as in our previous work, including low and 
intermediate (z < 1.5) redshift galaxies and M stars. Those 
objects resulting in spectra with too low S/N to be identi- 
fied as either LBGs or lower redshift objects were classified 
as 'unknown' (see Table |3]). An example of a LBG with a 
Grade A redshift from clear Lyman alpha emission is shown 
in Fig. [1] Results of our spectroscopy generally confirm our 
priority ranking and show contaminant discrimination is in- 
creasingly difficult from photometry alone for fainter ob- 
jects. Of the highest ranking objects, ~20 per cent are spec- 
troscopically confirmed as LBGs (see Table [3]). A significant 
fraction of the other priority 1 and 2 sources have no clear 
spectroscopic classification. In the following we discuss our 
results field-by-field. 



3.1 J0338+0021 

Our observations confirm six LBGs (four grade A, two grade 
B) in the field of J0338-f0021. Five of these cluster within 
Az < 0.05 (< 2500 kms'^) of the QSO redshift, with a sixth 
at Az = 0.06. From the ERGS survey, we would expect typ- 
ically zero or one object in a Az = 0.1 bin at this redshift 
in an unclustered field, given our observational setup. The 
LBGs do not spatially cluster around either QSO (Fig. [2}. In 
this field we also targeted a lower-priority object which was 
situated within a few arcsec of the QSO. Emission from the 
QSO increased the uncertainty on its photometry, contribut- 
ing to it being formally assigned a lower priority. However, 
its proximity to the QSO made it an intriguing candidate for 
spectroscopic follow up. Spectroscopy (Fig.[3)l showed a sin- 
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Table 1. The exposure times, seeing and 2a magnitude limits for the three QSO fields. 

Field Filter Total Exposure /hrs Average Seeing /" 2ct Magnitude Limits 

27.5 
26.9 
25.9 
27.6 
27.f 
26.2 
27.8 
27.f 
25.8 



J0338+0021 


R Special 


1.61 (20x290s) 


0.76 


" 


I Besael 


1.60 (24x240s) 


0.71 


II 


z Gunn 


2.17 (65xl20s) 


0.65 


Jf204-0021 


R Special 


2.09 (26x290s) 


0.85 


" 


I Bessel 


1.60 (24x240s) 


0.75 


" 


z Gunn 


2.00 (60xl20s) 


0.67 


J2f30+0026 


R Special 


1.53 (19x290s) 


0.78 


" 


I Bessel 


1.60 (24x240s) 


0.85 


" 


z Gunn 


2.33 (70xl20s) 


0.65 



Table 2. The number of priority 1 and priority 2 objects in each field and the fraction of these observed. Overall a high proportion of 
the total number of priority 1 and priority 2 targets were observed. 



Field 



J0338 
J1204 
J2130 



All 



Priority 1 
targets 



Priority 1 
on mask 



Priority 2 
targets 



42 

24 
27 



93 



20 (48%) 
18 (75%) 
12 (44%) 



50 (54%) 



32 
26 
46 



104 



Priority 2 
on mask 



13 (41%) 
15 (58%) 
15 (33%) 



43 (41%) 



gle, spatially-unresolved line at 7327 A coincident with the 
source position. The line is consistent with Lya at ^ = 5.024, 
within Az = 0.003 (Ai;=150 kms"^) of the QSO. This emis- 
sion line is highly unlikely to be due to extended Lya emis- 
sion from the QSO rather than Lya emission from an in- 
dependent source - firstly, because of the redshift offset be- 
tween th^_tw2_souJcesand_second^ experi- 
ence (jBremer et al.lll992l : iHeckman et aLlllOQin shows that 
extended QSO emission should be spatially resolved and fill 
at least part of the slit. 

Within this field we find nine objects, including the two 
QSOs, that have redshifts separated by a few thousand km 
s~^. Although the system is not likely to be viralised and 
the uncertainties are large given the small number of objects, 
the line-of-sight dispersion of the system can be estimated. 
We fit a Gaussian plus a constant model to the distribution 
of Lya determined redshifts, in order to model the galaxies 
in the redshift spike and any foreground/background galax- 
ies together. In this case the distribution is well fit by just 
a Gaussian, with almost no constant component to repre- 
sent foreground/background galaxies, giving a line-of-sight 
dispersion of the galaxies in the redshift spike (the Gaussian 
component) of 90014|!|q° kms~^. The errors are the la errors 
of the fitting procedure, which dominate over the redshift 
errors. 



3.2 J1204-0021 

Spectroscopy of the J1204-0021 field confirmed redshifts for 
four LBGs (three grade A, one grade B) from the high prior- 
ity targets over the two masks. Of the four spectroscopically- 
confirmed LBGs, two faU within Az < 0.04 (< 2000 kms"^) 
of the QSO redshift, a third at Az = 0.11 and the fourth 
clearly unconnected with the QSO at Az — 0.34. The two 
LBGs with redshifts close to that of the QSO are within 1 
arcminute of it on the sky (Fig. Q. Again, from the ERGS 
survey, we would expect zero or one object in a Az — 0.1 red- 
shift bin for a field showing no sign of clustering. Thus while 
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Figure 2. The 2D spatial distribution of LBGs relative to the 
J0338-I-0021 QSO (star). There is also a second, less luminous 
QSO near the edge of the field (star). The LBGs within Az = 
0.05 of the SDSS QSO are shown with squares and the other 
priority 1 targets from photometry are shown in the background 
as crosses. 



the significance of any clustering around the QSO in this 
field is obviously less than that for the J0338 field, confirm- 
ing two objects close to the QSO (in combined spatial pro- 
jection and redshift) implies some level of clustering around 
the QSO - a level that certainly could not be discerned from 
the numbers of LBG candidates selected from the photom- 
etry alone. 



3.3 J2130+0026 

Spectroscopy of the J2130-I-0026 field confirmed just two 
LBGs across the two masks. Neither were at redshifts close 
to that of the QSO (Fig. [5} . The imaging of this field con- 
tained significant scattered light from nearby bright stars, 
the effects of which were only properly understood and dealt 
with after the spectroscopic masks had been designed. Once 
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Table 3. The number of confirmed z^^S LBGs, low redshift objects and unknown sources due to low signal-to-noise spectra. Low-redshift 
objects have either clear continuum across the whole spectrum, an undulating spectrum indicative of a M-star or two or more emission 
lines indicating a low-redshift galaxy. Objects whose spectrum did not show evidence of being at low or high redshift were classed as 
unknown. Generally the higher priority objects were found to be LBGs and the lower priority objects to be low-redshift objects but there 
are also many priority 1 objects not confirmed as z~5 LEG, due to the low signal-to-noise spectra of faint sources making it difficult to 
see the continuum or be conclusive about a possible break. The higher priority objects are generally fainter, as they were selected to be 
at high redshift and hence their spectroscopy is more likely to have low signal-to-noise. 



Sample Confirmed high z 



Low z 



Unknown (low S/N) 



Total 



Pri 1 


6 {+2 QSOs) (18%) 


16 (36%) 


21 (46%) 


43 {+2 QSOs) 


Pri 2 


4 {+1 QSO) (7%) 


24 (42%) 


29 (51%) 


57 (+1 QSO) 


Pri 3 


1 (4%) 


19 (68%) 


8 (28%) 


28 


Pri 4 


(0%) 


34 (79%) 


9 (21%) 


43 



Table 4. The spectroscopically confirmed z~5 LBGs and QSOs in the three fields. The error on the redshift is 0.001 for grade A LBG s 
and 0.003 for grade B LBGs. The redshifts of the J0338+0021, J1204-0021 and J2130-I-0026 QSOs from SDSS (Schn eider et alj|2010t) . 
which uses absorption lines to determined the redshifts, are 5.0319, 5.0319 and 4.951 respectively. In all cases we have determined the 
redshifts from the peak of the Lyo line and hence the slight discrepancy between the observed redshift and the published redshift for 
the two QSOs we targetted. *Spectroscopy of these objects was not obtained and the redshifts stated below are from the literature 
toioreovski et alj|2003l : ISchneider et alj|2010h . tThe close companion to QSO J0338 (see text). Its photometry may be contaminated by 
scattered light from the QSO. 



Slit Rcf. 



RA /deg Dec /deg Redshift 



I /mag 



R-I /mag I-z /mag Priority Grade 



J0338-I-0021 QSO 


54.62211 


0.36553 


5.027 


19.9±0.001 


1.3±0.002 


0.1±0.003 


1 


A 


RD567 QSO* 


54.62513 


0.31122 


4.960 


21.1±0.003 


1.2±0.008 


O.OiO.Ol 


na 


na 


lall 




54.57256 


0.40150 


5.00 


24.9±0.13 


1.4±0.28 


0.1±0.33 


1 


B 


lal2 




54.58068 


0.40493 


5.090 


25.5±0.19 


1.7±0.42 


0.2±0.48 


1 


A 


2alt 




54.62303 


0.36554 


5.024 


25.7±0.46 


1.2±0.14 


0.1±0.30 


2 


A 


2a6 




54.65893 


0.38288 


5.027 


24.8±0.13 


2.0±0.39 


0.0±0.33 


1 


A 


2all 




54.64906 


0.40444 


4.987 


24.7±0.13 


2.1±0.39 


0.4±0.24 


2 


A 


2b3 




54.62834 


0.32185 


5.050 


23.7±0.03 


1.2±0.08 


0.1±0.09 


3 


B 


2bl0 




54.68122 


0.35007 


5.069 


24.9±0.13 


1.4±0.28 


0.9±0.18 


2 


A 


J1204-0021 


QSO 


181.17406 


-0.36364 


5.086 


19.2±0.0006 


1.5±0.002 


-O.litO.OOl 


1 


A 


IblO 




181.17325 


-0.38480 


5.424 


26.1±0.31 


>1.2±0.59 


<0.6±0.55 


1 


A 


lbl4 




181.18966 


-0.37056 


4.977 


24.8±0.13 


1.4±0.20 


-0.1±0.24 


1 


A 


2al 




181.18042 


-0.36293 


5.101 


25.5±0.20 


>1.8±0.54 


0.0±0.50 


1 


A 


2a3 




181.18167 


-0.35564 


5.050 


24.9±0.13 


1.5±0.26 


0.7±0.18 


2 


B 



J2130+0026 QSO* 
la5 
laS 



322.53727 0.43612 4.951 20.5±0.004 0.9±0.005 -0.3±0.009 na na 

322.50040 0.45221 4.777 25.2±0.15 1.0±0.22 -0.4±0.45 2 A 

322.50383 0.46411 5.133 25.4±0.23 >1.8±0.44 0.0±0.44 1 A 



0.30 ,"" * ' 
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Figure 3. The 2D and ID spectra of the close companion to the 
J0338 QSO. The 2D spectrum is smoothed by the 2 by 2 pixel 
boxcar average and the possible Lyman alpha emission is circled. 
Again the sky lines are plotted underneath for comparison. At an 
I magnitude of ~26, this source is too faint for any continuum to 
be detected. 



Figure 4. The 2D spatial distribution of LBGs (diamonds) rela- 
tive to the J1204-0021 QSO (star). The LBGs within Az = 0.05 
of the QSO are shown with squares and those further away by 
diamonds. The priority 1 targets from photometry are shown in 
the background as crosses. 
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Figure 5. The 2D spatial distribution of LBGs (diamonds) rela- 
tive to the J2130+0026 QSO (star). The LBGs within Az = 0.05 
of the QSO are shown with squares and those further away by 
diamonds. The priority 1 targets from photometry are shown in 
the background as crosses. 



Figure 6. The total number of LBGs found per 0.05 redshift bin 
away from the QSO. There is a clear spike at the redshift of the 
QSO but the majority of that signal comes from just one field, 
J0338. J2130 has no confirmed LBGs at the redshift of the QSO. 



the effects of the scattered light were fully corrected for, the 
subsequent photometry indicated that up to eight extra high 
priority objects could have been placed on the mask. With a 
spectroscopic confirmation rate of around 20% for high red- 
shift LBGs, this suggests up to 2 more LBGs may have been 
detected in this field. Therefore, it is possible that clustering 
similar to that seen in the J1204 field may be present here. 
Nevertheless, it is obvious that this field does not show the 
clear and strong clustering signature of the J0338 field. 



4 DISCUSSION 

4.1 Comparison to clustering in ERGS 

If we consider the spectroscopic results of the three fields to- 
gether (see Fig. |6]), it is clear that LBGs do cluster around 
2 ~ 5 SDSS QSOs. However, the strength of the clustering 
signal varies from QSO to QSO. Most of the signal in Fig. [6] 
comes from J0338 and none from J2130. While this is consis- 
tent with the often contradictory results found in the liter- 
atur e for the properties of high redshift QSO environments 



fe.g.lStiavelh et al.ll2005l : IZheng et al.]|2006l : iKim et al.ll2009l : 
lUtsumi et al.ll20ld ). the use of spectroscopy, along with the 
availability of an effective control sample from the ERGS 
survey, confirms that this variation is real rather than an ef- 
fect of different observational strategies and QSO redshifts. 
The comparison we have made takes into account both the 
completeness of the spectroscopic follow-up and the effec- 
tiveness of the spectroscopy, given that not every LBG can 
be spectroscopically confirmed with the quality of the data 
and the exposure times used. 

In particular, we can compare the clustering of LBGs 
around QSOs to that found in the ERGS fields - which 
were not selected on the presence of a known distant ob- 
ject. Eight of the ten fields in ERGS show no sign of clus- 
tering in redshift and, at 2 ~ 5, show typically zero or one 
spectroscopically-confirmed source per Az = 0.05. The same 
is true of the fiel d targeted with the s ame observational setup 
and strategy in iLehnert fc Bremer (12003') and three other 
fields that surround it ( Daviesll201ll ) . Additionally the three 



QSO fields in this work show no sign of any clustering at 
redshifts other than that of the QSOs. 

However, two fields in ERGS show significant spikes in 
their redshift distribution, one with seven objects between 
5.11 < z < 5.21 (J1040; see Fig. [8]) and another with 17 
objects between 4.95 < z < 5.15 - most of these split into 
two narrow spikes in the redshift distribution (J1054-12; see 
Fig- El- Plotting the spatial position of the objects in the two 
spikes in the J1054-12 field (Fig. llOp shows that they cluster 
spatially as well as in redshift. It appears that there are two 
separate structures in this field that may potentially merge 
to form a single larger one. If a representative a in km s~^ is 
determined for each of the ERGS structures using the same 
method as for the spike in the J0338 field, similar values are 
obtained: ~ 1300tlfo° km s"^ for J1040 and 850^3^^° and 
640l2?o° km s"^ for the two peaks in the J1054-12 field. 
This indicates that the redshift clustering in these fields is 
comparable or stronger to that in the J0338 field - although 
we note that if the redshift spikes in J1054-12 were not able 
to be separated in redshift, then we would have a broader 
distribution like that in J0338. However, the richest of the 
ERGs fields, J1054-12, was foUowed up with five MOS masks 
allowing a highly complete survey of all high priority objects. 
As with our other QSO fields, the J0338 field was observed 
with just two masks, in this case allowing about half of the 
highest priority objects to be observed (see Table [2]). It is 
likely that if more masks were used the number of objects 
in the redshift spike would increase, potentially to the level 
of that seen in J1054-12. 

Given these results, searching for high-redshift struc- 
tures by targeting QSOs appears to be not much more ef- 
ficient than observing blank fields. However, if we consider 
the statistics not in terms of fields, but in terms of inde- 
pendent redshift bins, the picture is subtly different. If we 
take all the 'blank' sky fields, from the ERGS survey and 
the fo ur fields in the area targeted by iLehnert fc Bremeij 
12003), ^iid randomly sample bins of width Az = 0.05 be- 
tween 4.8< z <5.3, we find only 7 per cent of such bins 
contain 2 or more LBGs. Similarly, 3 per cent of such bins 
contain 3 or more LBGs and 2 per cent contain 5 or more 
(see Fig. [7)l, with the majority of these bins arising from 
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Figure 7. The cumulative probability distribution of finding n or 
more (with n running from zero to eight) LBGs in a Az = 0.05 
redshift bin in the 14 FORS2 fields comprising the ERGS and 
BDF surveys. Most of the signal in the n ^ 2 bins arises from the 
two structures in Fig. [8] and [9] 



Figure 9. A histogram of the distribution of spectroscopically 
confirmed LBGs found in the J1054-12 field of ERGS. Shown 
underneath is the spectrum of sky emission. 



4.85<z<4.9 
4-95<z<5-0A 




Figure 8. A histogram of the distribution of spectroscopically 
confirmed LBGs found in the J1040 field of ERGS. Shown under- 
neath is the spectrum of sky emission. 



the two ERGS fields showing spikes in their redshift distri- 
bution. Given the results presented here, where two of the 
three QSOs show evidence of objects at the same redshift, 
LBGs tend to cluster around 2 ~ 5 QSOs more often than 
they do around an arbitrary point at high redshift. 

It appears that, as well as there being a clear varia- 
tion in QSO to QSO environment marked by surrounding 
LBGs, even the strongest clustering around z ^ 5 QSOs is 
no stronger than that seen in the field. Of the three QSO 
fields, one shows no evidence for clustering, another shows 
significant evidence for clustering of LBGs, whilst the third 
shows some evidence of clustering, but with only two LBGs 
at the QSO redshift. Given that the probability of finding 
two LBGs (three if the QSO host is included) in a single 
random redshift bin not associated with a richer structure 
is so low in our previous work (Fig. [7]), even two LBGs at 
the same redshift as a QSO indicates that the QSO is found 
in an over-dense structure. 



A 
A + 



163,6^ 165.65 163.58 153.70 163.72 163.7^ 



Figure 10. The 2D spatial distribution of spectroscopically con- 
firmed LBGs in the ERGS field J1054-12. Those in the two red- 
shift spikes are highlighted as diamonds and triangles and the 
dotted lines shows the extent of the FORS2 field of view. 



4.2 Clustering scales and markers 

Analysis of the Millennium-XXL simulation bv lAngulo et al.l 
(2012) showed that the eventual fate of halos hosting high- 
redshift QSOs (specifically those at z — 6) is most reliably 
determined from the richness of their environments on scales 
of 7-14 comoving Mpc. Those with the richest environments 
on these scales have the highest probability of ending up in 
massive clusters at the present epoch. Consequently, cluster- 
ing on scales of a few Mpc may not be the strongest signature 
of a forming galaxy cluster. 

Here we note that the over-density of LBGs in the J0338 
field is spread over (and probably extends beyond) the ~ 
7x7 arcmin^ (~ 15 x 15 comoving Mpc^) field, rather than 
clustering on smaller scales around just one of the QSOs. 
The over-densities in the ERGS fields are also spread over 
similarly lar ge scales, as is the lik ely z ~ 6 protocluster 
identified bv lToshikawa et al.l (|2012i '). rather than clustering 
on the 1-2 arcminute (~ 2 comoving Mpc) scales typical of 
more evolved clusters seen at z ^ 1 — 2 (3 Gyr later). Given 
the results of lAngulo et al.l l|2012r ) , it is likely that these 
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over-densities mark out the first stages in the formation of 
massive clusters. 

In our current analysis we have assumed that LBGs are 
reliable tracers of the underlying density distribution, and a 
lack of LBGs clustering at the redshift of the QSOs indicates 
that the QSO environment is not particularly over-dense. It 
is of course possible that we are being misled and that the 
number of detectable LBGs does not strongly correlate with 
the density field around a QSO. Given that galaxies selected 
through the Lyman break technique are simply those with 
significant unobscured rest-frame UV emission, arising from 
ongoing star formation (at rates of ~ 10 — 30 M0yr~^), it is 
unlikely that these are anything other than typical galaxies 
at their redshift. They should therefore trace the underly- 
ing density field as well as any other population which is 
undetected in the rest-frame UV to our flux limits. 

It is possible that the presence of the QSOs may in- 
fluence the colours and morphologies of galaxies in their 
immediate neighbourhoods due to an interaction between 
those galaxies and both the ionizing radiation fleld and mat- 
ter outflows from the AGN. Such galaxies would lie only a 
few arcseconds from the QSOs as we expect the AGN to 
be beamed in our direction, and for a galaxy to be caught 
within the beam it would need to be projected close to the 
QSO line of sight. Also it would need to be physically close 
to the QSO simply due to the drop-off with distance of the 
radiation fleld and outflow. We paid careful attention to all 
objects above our 7— band magnitude cut in close proximity 
to the QSOs, both because of the possibility of interaction 
and because the proximity to the QSO could increase the un- 
certainty on the photometry. Where slit contention allowed 
we placed them on the masks e.g. the close companion object 
to QSO J0338). Unless there is a population that is signifl- 
cantly reddened through interaction with a QSO (and there- 
fore drops below our magnitude cut) it is unlikely that such 
objects were preferentially excluded from our spectroscopy. 

Limited work has been carried out to explore the envi- 
ronments of high-redshift QSOs at other wavelengths, which 
might uncover redder or m ore obscured clustering popula- 
tions. IPriddev et al.l (12003) found an excess of S'ssom'" > 4 
mjy sources within ~ 1 arcminute of three z > 5 QSOs, in- 
dicating that dusty starbursts may cluster in their environ- 
ments (two of the QSOs themselves were detected as hyper- 
luminous dusty starb ursts). One of the Pr iddey QSOs was 
previously studied bv IStiavelli et al.l (|2005l ). who identifled 
tentative photometric evidence for an excess of distant LBGs 
in its immediate fl eld. As yet non e of th e sour ces in the three 
fields studied by iPriddev et al] l|2008f) and IStiavelli et al] 
((20051) have published spectroscopy confirming the excess 
objects to be at the redshift of their respective QSO. 

4.3 Is J0338 typical of z ~ 5 SDSS QSOs? 

As most of the clustering signal seen in this work arises 
from a single QSO field, it is worth asking whether or not 
that QSO is typical of the z ~ 5 SDSS QSOs. Ahhough 
not taken into a c count when this QSO field was selected, 
iDiorgovski et al.l (|2003l ') discovered a second fainter QSO 
(RD657) ~ 3' away from SDSS J0338-I-0021 with a very sim- 
ilar redshift {z = 5.02 for J0338, z = 4.96 for RD657). They 
argued that this pair marked a possible large-scale structure 
or protocluster. While there are no reported similar near- 



neighbours of the other two QSO fields probed in this work, 
there is no evidence in the literat ure that such neighbour s 
have been searched for. However, IDiorgovski et al.l (|2003l ) 
note that of the 13 other z > 4.8 QSOs studied by them, 
none appear to have a similar near-neighbour QSO. There- 
fore, it is reasonable to assume that the occurrence of QSO 
pairs on this scale is infrequent. Given the rarity of the QSO 
phenomenon, a rich environment clearly raises the probabil- 
ity of finding two QSOs in close proximity. 

QSO J0338-I-0021 is also a strong mm/submm source 
with continuum detections from 350^m to 1.2 mm, and 
detections of CO (5-4), [N il] and [CII] line emission (sum- 
marised in lWu et al.ll200 i/). indicating that the host galaxy 
is undergoing a strong (~ 2500Mp)y r~^) starbu rst. Of the 
other two QSOs in this work. ICarilli et al.l (|2001f ) measured 
the fiux of J1204-0021 to be S = 0.6 ± 0.4 mJy at 250GHz 
(1.2 mm), six times fainter than J0338+0021, and no mea- 
surement of J2130-t-0026 at mm/submm wavelengths has yet 
been published. 

A subset of other high redshift QSOs have been ob- 
served at mm and subm m wavelengths (e.g. IWillott et al.l 
l2007l : IPriddev et aIll2008lV Of those high redshift QSOs ob- 
served at 1.2 mm (e.g . ICarilli et al.l 1200 ll : iBertoldi et al.l 
l2003l : iPetric et al.ll2003l i. J0338-F0021 appears to be one of 
the most luminous (within the top 15 per cent) indicating it 
is among a subset that are hyperluminous, with dust masses 
of ~ 10* M© or more and star formation rates of > 1000 
M0 yr~^. The presence of a starburst of this nature is a 
strong indication that the QSO host is particularly mas- 
sive, and the more massive an early galaxy, the more likely 
it is to be found in a dense environment at high redshift 
(jSpringel et al.ll2006l ). Given the high SFR (and hence in- 
ferred mass) of this host galaxy, the high sub- mm brightness 
of the source may be a more important infiuence over the 
strength of any surrounding over-density than the fact that 
it is a QSO. 

Little observational evidence for the nature of the envi- 
ronme nts of very d is tant hyperluminous sources currently 
exists. ICapak et al.l (|201ll ) identified a likely protocluster 
containing both a QSO and a strong sub-mm luminous 
starburst at z ~ 5.3, whilst observations of the two most 
over-dense ERGS fields at mm and sub-mm wavelengths 
resulted in no detections, alth ough the full spatial extent 
of both fields was not p robed (|Stanwav et al.l l2008al . I2OI0I : 
iDavies et al.ll2010l . l2012l ). However, our observations in the 
two ERGS fields are sensitive enough to detect any submm 
sources as bright as QSO J0338+0021 o r indeed as br i ght a s 
~ 40 per cent of high redshift QSOs in lCarilh et"ai1 (1200 ll ) 
etc., suggesting the presence of a sub-mm source is not es- 
sential for an over-density of this richness to form. 

Inherent in the modelling of lAngulo et al.l (|2012l ) and 
other similar studies is an uncertainty in the fraction of mas- 
sive high redshift structures that contain luminous QSOs. In 
the above we have assumed that the two structures identi- 
fied in the ERGS fields differ from those around the QSOs 
in that they do not contain a powerful opticalUy-luminous 
AGN. However, it is possible that they host QSOs that are 
beamed away from us, QSOs that are currently in their "off' 
state or are otherwise obscured. If so, such sources may still 
be sub-mm bright objects even if the QSO emission is unde- 
tectable in the optical. Given our existing mm and sub-mm 
data for these fields, we have no evidence that the two struc- 
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tures marked out by LBGs in the ERGS fields are centred 
on a misdirected or "off' QSO. 



5 CONCLUSIONS 

We have carried out an imaging and spectroscopic survey of 
the near (~ 6 arcminute) environment of three z ~ 5 SDSS 
QSOs using a similar observational setup and selection to 
that used in ERGS. The similarity of approach confirmed 
that our strategy for finding LBGs is effective and provided 
examples of strong redshift clustering with which to compare 
the results in the QSO fields. 

We spectroscopically confirmed the presence of two or 
more LBGs at the same redshift as the QSOs in two of the 
three fields. Given that the probability of finding two or 
more LBGs in a narrow {Az ~ 0.05) redshift range in 14 
previously observed blank fields in ERGS and BDF is below 
7 per cent, it is clear that QSOs do reside in over-dense 
regions relative to the field at z ~ 5. It is also evident that 
there is significant variation in the strength of clustering of 
LBGs from QSO to QSO. Previous studies that relied on 
photometrically identifying over-densities of LBGs around 
QSOs would have easily missed the more weakly clustered 
fields, where spectroscopy is needed to show that the QSO 
and LBGs are found in the same structure. Therefore, null 
results in these previous studies may have underestimated 
the richness of the QSO environments. 

Even the richest of the three QSO environments studied 
is no richer than the discrete structures identified in blank 
sky fields of ERGS. When all blank sky fields observed by 
us are taken into account, one such structure is detected per 
~ 7 fields. If we split these fields into redshift bins of width 
Az = 0.05, only ~ 2 per cent contain structures as rich 
as that found in the richest QSO field. Hence, significant 
structure is more often found around QSOs than in 'blank' 
sky regions. 

We have been hesitant to call both the ERGS and the 
richest QSO-marked structures proto-clusters, as to do so 
requires evidence that the clustered regions will eventually 
form massive (M > lO^'^M©) structures at lower redshift. 
However, by the standard of other claims for the discovery 
of such systems that exist in the current literature, these 
structures, with velocity dispersions of < 2000 km s~^, are 



very likely to be 'proto-clusters' 
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